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1 Introduction 

The high precision of experimental measurements of benchmark processes such as jet, 
heavy quark, vector boson and Higgs boson production at the LHC calls for an equally 
precise theoretical evaluation of relevant production rates. In particular, exact perturbative 
predictions beyond next-to-leading order accuracy are desirable and sometimes essential 
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to reduce the theoretical uncertainty. Accordingly, in recent years considerable effort has 
beed devoted to computing next-to-next-to-leading order (NNLO) corrections to various 
production rates. Fully differential cross sections have been evaluated for vector boson 
[1, 2], Higgs boson [3, 4], diphoton [5] and Higgs-vector boson [6] production, and the 
computation of the total cross section for top-antitop production at NNLO is well underway 
[7-9]. Jet rates and event shapes in electron-positron annihilation to 2 and 3 jets have also 
been computed [10-18]. 

The construction of a general subtraction scheme to calculate cross sections at NNLO 
accuracy has also received significant attention [19-37]. Recall that the the NNLO correc- 
tion to the cross section for producing m jets is a sum of three contributions, the doubly 
real, real-virtual and doubly virtual terms, 



a 



NNLO 



/ da^ 2 J m+2 + I da m v +1 J m+ i + / d<x^ v J m , (1.1) 

J m+2 J m+1 Jm 



where each contribution on the right-hand side is separately infrared divergent in four 
space-time dimensions. The essential idea of subtraction is that by subtracting and adding 
back suitable approximate cross sections, we can reshuffle divergent pieces among the three 
terms in Eq. (1.1). This reshuffling amounts to writing Eq. (1.1) as follows, 
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NNLO 



m+2 



.NNLO , 
dcr m+2 + 



m+1 



j NNLO , 
dcr m+1 + 



da. 



NNLO 



where the three terms 
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da. 
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+ 



1 



j RV,Ai , 
dcJ m+l + 



, RR,A, 
dcJ m+2 



}. 



(1.2) 

(1.3) 
(1.4) 

Jm > 
\l.5) 

are each integrable in four dimensions by construction [20, 22, 23, 38]. Importantly, just 
as at NLO accuracy [27, 38-41], all approximate cross sections appearing in Eqs. (1.3)- 
(1.5) can be constructed once and for all, independently of the process or observable being 
studied. 

The approximate cross sections that contribute to dcr^S and d<T m+i° m -^ c l s - (1-3) 
and (1.4) were precisely defined in Refs. [22] and [23] respectively. The integral of the 
real-virtual approximate cross sections, denoted formally by L in Eq. (1.5), were computed 
in Refs. [25, 26, 28], while J 2 da^! 2 12 was evaluated in Ref. [32]. Finally, in Ref. [42], we 
calculated all collinear-type contributions to the integrated doubly unresolved approximate 
cross section, j 2 d<r^'^ 2 . 

In this paper, we complete the definition of the subtraction scheme by computing 
the double soft-type contributions to / 2 da m+2 2 - We use the method of Mellin-Barnes 
representations — developed in this context in Ref. [28] — to compute the various integrals 
that appear as building blocks for the complete expression. With the results of the present 
paper, the computation of the finite cross section in Eq. (1.5) becomes feasible. 
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2 Notation 



2.1 Matrix elements and cross sections 

Our notation was spelled out extensively previously [32, 42] and is recalled here only to 
the extent that we will need in this paper. 

The processes we consider are decays of a colourless initial state into any number of 
coloured patrons plus any number of additional non-coloured final state particles, that are 
however suppressed in the notation. We use letters form the middle of the alphabet, i, 
j, k, I, . . . , to denote resolved final state patrons, while r and s will label the two soft 
partons. The matrix element for a process involving n final state coloured partons will 
be denoted by \M n ), using the colour- and spin-state notation of Ref. [40]. The squared 
matrix element, 

\M n \ 2 = (M n \\M n ) , (2.1) 

is understood to be summed over colours and spins, which fixes the normalisation. In this 
paper, we will only need to use tree level matrix elements, which we denote |7Wi 0) ). 

Two- and four-parton colour correlated squared tree amplitudes are defined as follows 

\M^ {hk) \ 2 ^(M^\T t .T k \M^), (2.2) 

\ M n} im ,i)\ 2 = (M^m • T ki Tj ■ T t }\MV>) , (2.3) 

where Tj etc., is the usual colour-charge operator associated with the emission of a gluon 
from parton i. The square of the colour-charge operator, T 2 , depends only on the flavour 
of the emitting parton. We emphasise this by introducing the notation 

Cf^T 2 , (2.4) 

i.e., Cfc is the quadratic Casimir operator in the representation of parton i. Hence, C q = 
Cf = 1r(N^ — l)/N c and C g = Ca = 2Tr./V c . We also use squared colour-charge operators 
with multiple indices, e.g., T 2 r = Cf ir and T 2 rs = Cf irs . In such cases, the multiple index 
denotes a single parton with flavour obtained according to usual flavour summation rules: 
an odd (even) number of quarks and any number of gluons gives a quark (gluon). 

Finally, d<Jn\{p}) denotes the fully differential cross section for producing n partons 
at tree level with momenta {p} = {pi, . . . ,p n }. Its precise definition reads 




where M collects all factors that are independent of QCD; the symbol S{n} denotes sum- 
mation over different subprocesses; 5| n } is the Bose symmetry factor for identical particles 
in the final state; and d(p n ({p}; Q) is the phase space associated with final state momenta 
{p} = {pi, . . . ,p n } and total momentum Q 

dMpu . . . ,p n ; q) = n j2^ 5+{p * ] {27r)d5id) [ Q ~J^ p j- (2 - 6) 
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In this paper, we will use the cross sections for producing m and m + 2 partons, which we 
call the Born and doubly real contributions, respectively 

daf n ({p}) = d<4 0) (M) and da™ 2 ({p}) = dcrffi. 2 ({p}) ■ (2-7) 
2.2 Phase-space factorisation and kinematic variables 

The double soft-type subtraction terms are all defined precisely via the double soft mo- 
mentum mapping of Ref. [22], which maps the original set of m + 2 momenta into a set of 
m tilded momenta 

{p}m+2 ^ iP}i S) ■ (2-8) 

The explicit form of this mapping can be found in Ref. [22]. Here we recall only that it 
leads to an exact factorisation of phase space in the following form, 

d4> m+2 ({p} m+2 ;Q) = d0 m ({p}^ s) ;Q)[dpS(Pr,P,;Q)] , (2.9) 
where the factorised phase-space measure is 

[d*G(pr,p.;Q)] = d(\l)(\l)^- 1 ^- 1 ^dUpr,Ps,K-,Q)e(x rs )e(i - x rs ) , (2.10) 

with a massive momentum such that K 2 = X 2 S Q 2 . 

We will have use for the iterated double soft momentum mapping as well, 

{pU2^{?£i^{^ s) > (2-11) 

also defined in Ref. [22]. This mapping again leads to the exact factorisation of phase space, 
now in the iterated form 

d0 m+2 (W m+2 ; Q) = d^ m ({p}g^;Q)[dpfl(p r ; Q^dpg^fe; Q)\ . (2.12) 

The factorised phase-space measures read 

[dpfliPf, Q)) = dy tQ (l - ytQ^^-^dMPt, K- Q)Q{y tQ )Q{\ - y tQ ) , (2.13) 

where t = s and n = m + lori = f and n = m. The two-parton phase space dfoiPt, K\Q) 
can be written as follows in the rest frame of Q, 

dMPt, K- Q) = 2^ 1+2e vrT(l - e) de t e\- 2e 8{e t - y tQ )dn d ^(t) , (2.14) 

where St = 2Etj sfQ 2 is the scaled energy of parton t in the Q rest frame, while dO^_i (t) 
is the angular measure in d — 1 dimensions in this frame. In Eq. (2.14), S e is 1 



1 S e = (47r) E e _E7E as usually defined in the literature, and it is this factor which is conventionally included 
in the definition of the running coupling in the MS renormalisation scheme. Both definitions lead to the same 
expressions in an NLO calculation. At NNLO, the two definitions lead to a slightly different bookkeeping 
of IR and UV poles in intermediate steps, but physical cross sections are the same. Our definition leads to 
slightly simpler expressions at NNLO. 
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Finally, we recall the definitions of kinematic variables. First, two-particle invariants, 
such as Sjfc, SiQ, SfQ, etc., always denote twice the dot product of two momenta. A double 
index enclosed in brackets, as in e.g., Sj( rs ), indicates a sum of momenta. Hence 

s ik = 2pi ■ p k , SiQ = 2pi-Q, s f .Q = 2p r -Q, s i(rs ) = 2p { ■ (p r + p s ) , etc. (2.16) 

When these invariants are normalised to the total incoming momentum squared, we use 
the notation = sik/Q 2 ■ The momentum fractions for double and triple parton splitting, 
Zi )V and Zi jrs , are defined as follows 

Zi,r = — and z iyTS = lQ , (2.17) 

SiQ + S r Q SiQ + S r Q + S S Q 

with z rt i as well z r ^ s and z s ^ T given by obvious cyclic permutations. Clearly Zi^ r + z r ^ = 1 
and Zi. rs + z r ^i s + z Sj i r = 1. Furthermore, Sik(r) denotes the eikonal factor, 

Slk {r) = ^. (2.18) 

Lastly, the integrated counterterms turn out to be functions of the following two types of 
variables 



2J hr- and Y-rr Q = 

Q 2 ,Q y~i Q y~kQ 



Y,, t) (2.19) 



3 The integrated doubly unresolved approximate cross section 

We recall from Ref. [42] that after summing over the flavours of unobserved partons, the 
integrated doubly real approximate cross section can be written as 



^da^=da*®I$>\{p} m ;e). 



(3.1) 



Here the insertion operator 1^ reads 



4 O) ({p}m;<0 



+EE 

I iH 



5 (o),(j,0 



C A +^CSg' W) T 2 



TjT t (3.2) 



+ ^$MWM{T i T k ,T j T l }[, 

i,j,k,l ) 
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with 



p(0) _ 
^2,i — 


/ P (0)\ 
^irs ^ 


p(0) _ 


/p(0) 
1 iry's 


Cs(o).C?'.0 = 




q(O),(j,0 _ 
s 2 — 


\ D rs J 


q(O),(<,fc)(j,0 _ 
S 2 — 


(q(0)) { 

\ D rs J 



fi 



p pq(0) 

Ji 

p pa(0) 



p n(0)\ 

' fifj 



cs iria s#) 



(3,0 



(3.3) 



/< 

(7.0 

(i,fe)(i,0 



The flavour-summed functions, ( C-^ J etc., appearing on the right-hand side of Eq. (3.3) 

^ ' fi 

are given in terms of non flavour-summed functions in Ref. [42], where we also computed 
all those functions in Eq. (3.3) which do not involve double soft contributions. 

In this paper, we consider the remaining flavour-summed double soft functions. These 
are obtained from the non flavour-summed integrated subtraction terms as follows [42]: 



1. Triple collinear - double soft: 



p q(0) 



fll 



2. Triple collinear - soft collinear - double soft: 

(c ira cs ir;a s(°;) f = [C ir ,CS 

i.e., it is independent of the flavour /. 

3. Double collinear - double soft: 

i.e., it is independent of the flavours f± and /2- 

4. Double soft: 

,(<,fc)C?,0 1 



ir;s^rs J ' 



(3.4) 



(3.5) 



(3.6) 



g(0) 



g(0) 



2 <-'~' rs ! 99 

m _ i 

~ 2 



CM) + „, ra(o)iO',0 

rs Igg ' '"I IVrs Iqq 



(3.7) 



5. Soft collinear - double soft: 



i.e., it is independent of the flavour / 



(7,0 



rpo 



0,0 



(3.8) 
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4 Integrated double soft-type counterterms 



In this section we define explicitly all functions appearing on the right-hand sides of 
Eqs. (3.4)~(3.8), and compute them in terms of a basic set of integrals. 

1. Double soft: 



(i.fc)0'.0 



1 



(i,k) 

rs Igg 



S t / J 2 

X f(V0, VrQ + VsQ ~ Vrs, d' (in, €)) , 

(47F) -Q 2 ^ / [dpg(p r ,p s ; Q)] \s ik (r, s) 



S e / J2 

X f(V0, VrQ + VsQ ~ Vrs, d' (jn, S)) , 



Sir^ks H~ Skr$is Sik^rs ^.^ir^is 



(4.1) 



(4.2) 



x 



'i(rs)^fe(rs) 



'i(rs) 



7(2/0, J/r-Q + y s Q - Vrs, d'(m, e)) . 

(4.3) 



The non-abelian part of the double soft gluon current, Si k (r,s) in Eq. (4.2), reads 



5 ifc (r, S ) = 5^(r, S ) + 4 



(s.o.) / „\ i a ^ir^ks Sis^kr 



^i(rs)$k(rs) 



1 - e 1 



Js.o.) 



(r,s) 



-S ik (rs) , (4.4) 



where 



s£°-\r, s) = S lk (s) (S ls (r) + S ks (r) ~ S ik (r)) 



(4.5) 



is the form of this function in the strongly-ordered limit and Si k (rs) is given by 



S ik (rs) 



2s 



if. 



(4.6) 



Finally, the function f(y , y rQ + y sQ - y rs , d'(m, e)) appearing in Eqs. (4.1)-(4.3) rep- 
resents a small but convenient modification of the subtraction terms as compared to 
their original definitions in Ref. [22]. Its precise role and form are explained in Ap- 
pendix A, and in the following, we will include this factor without further comment. 

In this paper, we do not discuss the case when i, k, j and I in Eq. (4.1) are all distinct, 
this requiring at least four jets at NNLO. For the specific cases of two and three hard 
partons, the integrated double soft counterterms are computed in Appendix B and 
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we find 



(i,k)(j,k) 
(i,k){i,k) 



r s (0)i(i,fe) 

L^-Ys Igg 



(i L k) _ Tr_ 
s m - Ca 



l2S,i(Y 1%Q , Y~ ilQ , Y 1%Q , e; y , d' ) , (4.7) 

l2S,2(Y 1%Q ,e;y ,d' ), (4.8) 

l2S,2(Y-r ^ >Q ,e;y ,d' Q ) -X 2 s,3(Yr^ Q ,e;yo,d' ) 

2X 25 ,4 (Yr £ Q , e; y , d' ) + l 2 s,5 (Yr ^ Q , e; y , d' ) 

4Z2Sfi(Yr% Q ,e;y ,d' ) + 2(1 - e^S/K^^g, e; 2/0,4) 

(1 - e)X 25 , 8 (e; y ,d' ), (4.9) 

■^s.e^fc^^^o^o) - 2X2 5 ,7('5 / 7fc j Q,e; yo,d' ) 

+ 2X25,9(6; y , 4)1 ■ (4-10) 



The integrals ^25^ for n = 1, . . . , 9 are denned in Eqs. (B.6)-(B.14). 



2. Triple collinear - double soft: We decompose the triple collinear - double soft gluon 
counterterm as a sum of abelian and non-abelian pieces, 



\c s(°)u =\c, q(°)l( ab ) + rn q(°)i( nab ) 

l KJ irs'- J rs ijgg l^ J trs'- J rs Igg ' l^ J trs'- J rs ifgg 



(4.11) 



following the decomposition of the triple Altarelli-Parisi splitting kernels in a similar 
fashion [43]. Then we have 



r C s (o)i(ab) = (i^) 2 ^ 



x f(yo, VrQ + y S Q - y rs , d'(m, e)) , 



4z? 



(4.12) 



\C,:„, 



(nab) 
irs^rs Ifgg 



[ d P { im(Pr,Ps;Q)} 



/ J 2 

(1 e) (Si r Z Sj i r SisZr,is) 



C A 
Tl 



ITS 



S; 



i(rs 



2zf 



2 1 



?.,rs 



+ 



+ 



+ (r <R- s) 



X /(Z/0, 2/rQ + - d' (m, e)) , 



(4.13) 



(4tt) 



-Q 



2< 



[ d i>£m (Pr , Ps! Q)] ^ S~ 
2 *j(rs) b rs 1 irs 



"i,rs 



Sis Zr 



•r,is 



~i~ %s,ir Si(rs)^rs{Zr,is ~i~ Zs,ir)^ 



x f(yo, VrQ + y S Q - y rs , d'(rn, e)) 



(4.14) 
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The integrated triple collinear - double soft counterterms are computed in Ap- 
pendix C. In terms of the basic integrals introduced in Eqs. (C.3)-(C.ll), we find 

ICirs^Jif = 4Z 2SAQ (e;y ,d' ) , (4.15) 
[C irs SW& b) = Sr- [ 4 ( X " ^25,11(6; yo, d' ) - 4(1 - e)X 25 , 12 (e; y , d' ) 

irs 

- 8X 25 , i3(e; yo, d' ) + 2X25, 14 (e; yo, d' Q ) - 4Z 2 5,i5(e; yo, d' ) 

- 2X 25 , ie(e; yo, d' ) + 2X 2(Sil7 (e; yo, d' ) + 2Z 2<5>:L8 (e; yo, d' ) , 

(4.16) 



[C irs Si°J] q q = —5- - 4Z 2 5,ii( e ;yo,d ) +4Z 2 5,i2(e;yo,c?o) + 2X 25 ,i3(e; yo, d' ) 



3. Double collinear - double soft 

^2 \ 2 



(4.17) 



(4.18) 



[<WW>] = (^Q 2 ') [ [dp%%r, Ps ; Q)} f ^ 

\ "e / J 2 °ir *T,i °js *aj 

x /(yo, 2/rQ + y S Q - y rs , d'(m, e)) . 
This integral is computed in Appendix D. The result is 

[C ir . js S(°J] = 4Z2 5 ,i 9 (Y TLQ ,e;y ,<,) , (4.19) 
where I 2 s,i9 is defined in Eq. (D.2). 

4. Soft collinear - double soft: 

x /(yo, yrQ + y S Q - y rs , d'(m, e)) . 

We remind the reader that above j / I, but (ir) may be equal to j or I. The cases 
when j, I 7^ (ir) and when e.g., j = (ir) lead to different integrals, as discussed in 
Appendix E. In terms of the basic integrals introduced in Eqs. (E.3) and (E.4), the 
result is 

[CS ir .. s S$] = -2X 2<s>20 (r IT)Q , e; y , d' ) , (4.21) 
if (ir) is distinct form both j and I, while for e.g., j = (ir), we find 

[CS ir;s S$]W) = -21 2S , 21 (xi,Y~ il>Q ,e;y ,d' ) . (4.22) 

5. Triple collinear - soft collinear - double soft: 

\2 



2 

r.is 



Sir ^(ir)s %r,is %s,ir (4.23) 

x /(yo, y r Q + y S Q - y rs , d'(m, e)) . 

This integral is computed in Appendix F. We find 

[C irs CS ir;s S( s )] = 4X25,22(^,6^0,4) + 4Z 25 ,23(x T ,e;y ,d , ) , (4.24) 
with X 2 s,22 and I 2 s,23 defined in Eqs. (F.3) and (F.4). 
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5 Computing the double soft integrals 



5.1 Simplifying the integrals 

Consider a generic double soft-type master integral, 

— <5 2e J J^P2^} F n{{sji,Zjy,e)f{yo,y rQ + y s Q-y rs ,d'(m,e)), (5.1) 

where .F n ({sjz, Zj}; e) is the integrand, which will depend on two-particle invariants Sji and 
possibly also momentum fractions Zj. (We use Zj generically to denote both the two- and 
three-parton momentum fractions defined in Eq. (2.17).) A straightforward treatment of 
Eq. (5.1) proceeds to first express the integrand in terms of independent (tilded) momenta. 
Then, choosing some particular Lorentz frame, the integral is written in terms of e.g., the 
angles and energies of momenta p^ and ps in this frame. The result obtained however 
turns out to be completely unwieldy. In particular, the extraction of e poles via sector 
decomposition is not possible, because one typically finds complicated singularities inside 
the domain of integration. It also turns out that such a parametrisation is not a useful 
starting point for deriving Mellin-Barnes representations. 

A much more manageable form is obtained, however, by observing that the value of 
the integral in Eq. (5.1) is unchanged if we replace the double soft unresolved phase-space 
measure [dp!^] by the iterated single soft phase-space measure, e.g., [dp^][dp^ +1 ]. 

To understand why this may be the case, note that by using the precise definitions of 
the single and double soft momentum mappings [22], it is easy to show that the two sets 
of tilded momenta {p}m^ and {p}m^ , where 

{p} m+ 2 ^ {P)t s) and {p} rn+2 A {p}$ +1 A , (5.2) 

are related simply by a three-dimensional rotation. This implies that the dot products of 
momenta in each set are the same for both sets. But the integrated counterterms depend 
only on these dot products (as opposed to the overall orientation of the tilded momenta), 
hence we are free to use the more convenient iterated phase-space mapping when computing 
the integrals, 

[dptm] F n{{sjh Zj}; e)f(y , y rQ + y sQ - y rs ,d'(m, e)) 

(5.3) 

( "\ ( \ 

5 - m +i\ Fn({sji, zj}; e)f(y ,y r Q + y s Q - y rs , d'(m, e)) . 

2 

By using the iterated form of the phase-space measure, the integrations over the variables 
of pr and ps may be performed sequentially. 

5.2 Computing the integrals via Mellin Barnes representations 

Following the discussion of the previous section, we write the generic master integral in 
Eq. (5.1) as 

l2S,n = {^Q 2e ) jf[d&][dpS^ (5-4) 
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Since the dimension of F n is [F n ] = (Q 2 ) 2 , we have 

F n ({sj h Zj}) = j^yF n ({y jU z jM }) , (5.5) 
and hence, using Eqs. (2.13) and (2.14), we find 

l2S,n = 2-^vr- 2 + 2 T 2 (l - e) £ dy rQ y^(l - y fQ ) d '^ f dy sQ y^(l - y sQ ) d '^ 
x ® (yo - VtQ - VsQ + yfQVsQ) dn d - 1 (r)dn d - 1 (s)F n ({y j i, Zj}) . 

(5.6) 

In the above form, Eq. (5.6) is directly suitable for treatment by Mellin-Barnes techniques. 
In particular, the angular integrations over the directions of pr and pt may be preformed 
sequentially using the results of Ref. [44]. Hence, it becomes essentially straightforward to 
derive a Mellin-Barnes representation for the integral: 2 

1. Substitute the specific form of F n . 

2. Perform the angular integration over dfirf_i(s), using the results of Ref. [44]. The 
momenta which are to be kept fixed during this integration are the intermediate, 
hatted momenta of Eq. (2.11). Thus, at this stage, all invariants need to be expressed 
in terms of these: 

Vik = 0- -y S Q)y~ik, Vis = yi s , y%Q = (i - ysQ)y\ Q + y~ is , i,k^s. (5.7) 

The result after this angular integration depends on (scaled) dot products of hatted 
momenta and perhaps Q, i.e., y* f , y~ iQ , etc. 

3. Repeat the previous step for the angular integrals over df^„i(f). The independent 
momenta with regard to the pr integration are now the final set of tilded momenta 
in Eq. (2.11), so we must express all invariants in terms of these: 

vik = i 1 - yfQ)nk > yir = yir> wq = i 1 - y?Q)y~ iQ + r ir , i,k^r. (5.8) 

After this integration, the expression depends only on invariants involving tilded 
momenta, as well as y r Q and y s Q, which are however integration variables themselves. 

4. Finally, perform the integrations over y r Q and y s Q to obtain the full Mellin-Barnes 
representation. If y$ = 1, the argument of the function in Eq. (5.6) factories and 
these last two integrations are automatically in the form of a beta function integral. 
For < yo < 1, additional Mellin-Barnes integrations must be introduced to bring 
them to this form. 

2 Note that contrary to the phase space integrals considered here, Mellin-Barnes representations for loop 
Feynman integrals can be constructed automatically, e.g., with the AMBRE.m package [45]. 
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With this procedure, we obtain a multi-dimensional Mellin-Barnes integral representa- 
tion of Eq. (5.6). It is usually not possible to evaluate these integrals analytically. However, 
for practical purposes, we are actually interested in the e-expansion of l2S,n, rather than 
the all-orders result. Importantly, this expansion may be performed in an algorithmic way, 
e.g., with the MB.m [46] or MBresolve .m [47] packages. Then we obtain the expansion 
coefficients as finite Mellin-Barnes integrals. We are able to evaluate all integrals that con- 
tribute to the 1/e 4 and 1/e 3 poles analytically, 3 while the rest of the expansion coefficients 
can be computed by direct numerical integration of the Mellin-Barnes representation. 

Finally, we call attention to the following technical detail. We find that it is best to fix 
the specific value of d' (see Eq. (A. 3)) before the e-expansion of the integrals. Therefore, 
in the following, we set d' Q = 3 — 3e, i.e. D' Q = 3 and d[ = —3. This is not an essential 
restriction and recomputing the expansions for different values of D' and d[ is in principle 
straightforward. 



6 Results 

6.1 Analytic expressions to 0(e" 2 ) 

We have obtained analytic expressions for the integrated counterterms up to 0(e~ 2 ) ac- 
curacy. The results below are quoted for d' = 3 — 3e, i.e., D' = 3 and d[ = —3 (see 
Eq. (A.3)), and generic yo £ (0, 1]. Before presenting the actual formulae, we call attention 
to the following features. Up to this pole order: 

• All kinematic dependence enters through logarithms of the variables x% and/or Y^q 
defined in Eq. (2.19). 

• Dependence on the cut parameter uq always enters in the same functional form, as 
follows, 

lny -2yo + y- (6.1) 

We observe that this is simply £(yo, 2) (recall that throughout we use D' = 3), with 
the E(z,N) function of Refs. [25] and [32], 

£(*, N) = In z - ~ k ~ ■ (6-2) 
k=i 

• Additional constants that enter are always recognised to be slight generalisations 
of jq/C q and J g /C g (since our integrated counterterms are dimensionless in colour 
space), where 7 9 = |Cf and 7 5 = ^-Ca — §Tkrif [48]. In order to emphasise this 
connection, we define the following two functions of nf [42], 

lq (nf) = \ and 7s (nf) = y - j^ n f • ( 6 - 3 ) 



3 In these calculations, the barnesroutines .m package of D. Kosower is useful for applying the first and 
second Barnes lemmas on Mellin-Barnes integrals automatically. 
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The formal 7if dependence of j q (rif ) is introduced in order to make possible a flavour- 
independent notation in the following. 

We find 

1. Triple collinear - double soft: 



-3-) ^-4s(y ,2)V^7 5 (n f ) 



2C /i Me 4 e 3 ' 



Si 



+ 0( e ~ 2 ). 
(6.4) 



2. Triple collinear - soft collinear - double soft: 

" 2 1 



C irs CS ir . s S^,J^j ^ (xi) — ^ 
3. Double collinear - double soft: 



e 4 £ 3 



lnx» + 7S(y ,2) 



+ 0(e" 2 ). (6.5) 



1 4 



S(y 0) 2) 



+ 0(6- 2 ). 



(6.6) 



4. Double soft: 



,(i,fc)0',fc) i 



^--3 ( lny ijfc)Q + ln^ feiQ + 4S(yo,2) 



2 



+ 0(e 

(note the lack of lij,Q dependence to this order) 

1 1 



(6.7) 



e 4 £ 3 



21ny ifciQ + 4S(y ,2) 



+ o( e ~ 2 ), 

(6.8) 



and 



S<$) (i ' k) (Y lk , Q ) = - 1 - 



+ 0(e" 2 ). (6.9) 



5. Soft collinear - double soft: 



^ ' Ji 



+ 0(e" 



(6.10) 



for j, Z ^ i, and 



7 ~ ^ 



+ 0(e" 2 ) 
(6.11) 



for e.g. j = i. 
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Substituting Eqs. (6.4)-(6.11) and the corresponding results from Ref. [42] — recalled 
here in Appendix G for the convenience of the reader — into Eq. (3.3), we obtain explicit 
expressions for the kinematics-dependent functions entering the insertion operator. We 
note that since in this paper we set D' Q = 3, for consistency we must also use the results of 
Ref. [42] with D' Q = 3. The results do not depend on do and «o U P to this order. Starting 



with eg , we find 



c$(*i) = -£sli + £r 

,(0) 



C A 

Cf, 



J 21nxi - 2S(y ,2) - 7/ . ("f^rj 



For C, we obtain 



21nxi -7/. (n f ) ) - hhiXj - j fj (n f ) 



+ 0(e" 



+ 0( £ " 2 ) 



(6.12) 



.13) 



However, eg- enters the insertion operator jg summed over its flavour indices, thus we 
are free to symmetrise in i and j. In particular, setting 

C^(xi, Xj ) = \(c^{x u x 3 ) + cW^Xij) , (6.14) 

(we use the usual notation of round brackets around indices to denote symmetrisation) we 
see that up to this order in the e expansion, cl°). simply vanishes, 



,(o) 



2,(«) 

(xi,xj) = 0(e" 2 ) 



The two-parton colour-correlated soft function, S 



(O),(i,0 



is 



I - i f 2 hi + 4S(y , 2) - 7s (nf ) J 



+ 0(6- 2 ) 



(6.15) 



(6.16) 



For CSg we obtain (note the vanishing of Yji,Q dependence to this order) 

CSg! ,(,V) (xi, = ^ (2 In - 2S(y , 2) - 7/j (n f ) ) + 0( e ~ 2 ) (6.17) 



when j and I are distinct form i, and 



CSg' (i,0 te,r«,Q) = ( 21nx, -2S(y ,2) - 7/i (n f ) ) + 0(e" 2 ) 



4e 3 



.18) 



for e.g. j = i. Last, we consider the four-parton colour-connected soft function, g(°)(*''W" ,i ). 
We find (note that Yi^Q dependence is absent to this order) 



3 (0),(i,fc)(j,fe) 



{Yik,Q,Yij,Q,Yjk,Q) 



1 1 



In Y iKQ + In y jfciQ + 4S(y ,2) 



+ o( e ~ 2 ), 

(6.19) 



if three indices are distinct, while 



$V m{l *\Y ik , Q ) 



1 1 



21ny ifciQ + 4S(y 0) 2) 



+ 0(e~ 2 ) 



.20) 



if only two indices are different. 
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6.2 Insertion operator for two- and three-jet production 



Using the results of Ref. [42] and the present paper, we can assemble the complete insertion 
operator relevant for two- and three-jet production. 

Let us consider first the process e + e~ — > 2 jets. The corresponding squared matrix 
element at tree level is \M^\l q , 2^) | 2 , i.e., the quark carries label 1 and the antiquark label 
2. Both the colour algebra and kinematics are trivial. Colour conservation implies 



TiT 2 



-C F 



(6.21) 



and Cf 1 = Cf 2 = C F , while momentum conservation gives 



Xl=X 2 = Y 12 ,Q = J/12 = 1 ■ 



(6.22) 



Then the insertion operator, Eq. (3.2), is a scalar in colour space and reads 



r(0) 



(^)] 2 [2C|(CS + C^-2CS^' 2) +4S 



(0),(1,2)(1,2) 
2 



2C F C A sj i 0) < (1 ' 2) 



5.23) 



with all arguments being equal to one. Using Eqs. (6.12)-(6.20), we find 



2v e Vq 2 



4 + x 72 + 29x - Ayrif 
~2^~ + 12p 



+ 01 



-2 



(6.24) 



where we used the notation [49] 



Ca 
C F 



Ik 
c F 



(6.25) 



Notice that up to this order, 1^ is independent of yo. 

The remaining expansion coefficients are computed numerically. We present these in 
Table 1 for the quantity 1% (pi,P2\e)/C%. In this calculation, we have set do = df = 3 — 3e 
and ao = yo = 1) however colour factors and the number of light flavours, rif, were kept 
symbolic. Hence the actual value of the expansion coefficient at a given order in e is given 
by the scalar product of the vector of numbers in the appropriate column of Table 1 with 
the vector of coefficients forming the first column. 

Turning to the process e + e~ — > 3 jets, the corresponding squared matrix element at 
tree level is \M.^ 2g, 3 9 )| 2 , i.e., the quark carries label 1, the antiquark label 2 and the 
gluon carries label 3. The colour algebra is still trivial and using colour conservation we 
find 

Ca - 2Cf m m mm Ca 

2 ' 



Ti To 



Ti T, = To T, 



(6.26) 
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coeff. 


1/6 2 


l/e 1 


finite 


1 


-7.416 ± 0.011 


-81.383 ± 0.052 


-236.572 ± 0.328 


C A 
c F 


-13.028 ± 0.003 


-66.274 ± 0.025 


-281.208 ± 0.188 


n f T R 
C F 


1.000 ± 0.001 


9.984 ± 0.010 


57.535 ± 0.075 



Table 1. Coefficients of the Laurent expansion 

of 2 V,2 for e+e ~ -> 2 J et production. We 

used do = d' — 3 — 3e and cto = yo = 1. 



Cj 2 = Cf, while C/ 3 = Ca- The insertion operator, Eq. (3.2), is again a scalar 



and Cf x 
in colour space and we have 



r(0) 



^S f ( ^ 



2vr 



Q 2 



n 2 



{^(cfjcxo + cgj^+^cgds) 



+ ^(4 ? (, 1) x 2 )+C^(x 2) x 1 ) 
+ CtC A (cgJ fl (xi,X3) + cg,(x 3 ,xi) 



(o) 



+ C 



(o) 



(x 2 ,X 3 ) + C^(X3,X 2 )) 



+ (C A - 2C F ) 



c A (csg^' 2 )(, 3 ,y 12 ) 
+ c F (cs^ (1 ' 2) (x 1 ,y 12 ) 



+ sf^(y 12 ; 



+ CE^ ( x 2 ,y 12 ; 



C A 



c A (c5^ (1 ' 3) (x3,y 13 ) 



+ s( 0) ' (1 ' 3) (r 13 ) 



cpfes^^^^^ + cs^^^,^) 



^ (0) ' (2 ' 3) (x 3 ,y 23 ) + sf'^(y 23 ; 



c F (csS'^(, 1 ,y 23 ) + csa^(, 2 ,y 2 3; 



2Ca(2C f - Ca) 



2 

f s 



,(0),(1,2)(1,3) 



(y 12 ,y 23 ,y 13 ) + s" 



(0),(1,2)(2,3) 



(yi2,y 3 ,y 



23 I 



2 

+ 2ci (sf ■(w- 3 ) (y 13) y 2 , r 23 ) + s?^ 8 ^ (y 23 , y 12j y 13 ; 

+ 2(2C F -C A ) 2 sf' (1 ' 2)(1 ' 2) (y i; 



(0),(2,3)(l,2)/ v y y V q (0),(l,3)(l J 2)/ V ; y y - 

1*23, *13, *12j + ^2 (*13, *23> *12, 



■12 



+ 2Ci(S 



(0),(l,3)(l,3)/ v s q (0),(2,3)(2,3) r y - 
2 1*13 J + >5 2 (*23, 



(6.27) 
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coeff. 


1/6 2 


l/e 1 


finite 


1 


17.072 ± 0.016 


-30.766 ± 0.088 


-253.694 ± 0.474 


C A 
c F 


36.428 ± 0.013 


35.832 ± 0.068 


-83.945 ± 0.294 


r<2 

C A 
Cp" 


12.580 ± 0.009 


3.999 ± 0.052 


-99.299 ± 0.296 


n f T R 

c F 


-16.945 ± 0.003 


-59.501 ± 0.023 


-128.771 ± 0.144 


C A n f T R 
r 1 


-3.879 ± 0.006 


-14.105 ± 0.052 


-29.974 ± 0.378 



Table 2. Coefficients of the Laurent expansion of — - — ' ^ 3 — for e + e — > 3 jet production. 
We used do = d' = 3 — 3e and ato = yo = 1. 



Substituting Eqs. (6.12)-(6.20) into this expression, we obtain 



2vr" 



Q 2 



eg 



n 5 3 2 

2 + -x+-x 2 
2 4 



+ 



77 2 

H a; 

24 



1 



-xyrif 



A + x 



-x 



6 + l2- X 



In yi2 



(6.28) 



3 2 
2x + -x^ 



(lnyi 3 + lny 2 3) 



where we used In I^q = In — In rrj — In . Note that the dependence on yo again cancels 
up to this order in the expansion. 

The rest of the coefficients in the e-expansion are computed numerically. By way of 
illustration, we present results for the fully symmetric configuration of final state momenta 
in Table 2. In this phase-space point, the various invariants take the following values 



2 

3 ' 



Yik,Q = 2 



and 



Vik — g 



(6.29) 



where i,k = 1,2,3 and i ^ k. We have again used do = d' Q = 3 — 3e and ao = yo = 1 during 
the calculation. As for the 2-jet case, colour factors and the number of light flavours were 
kept symbolic. The actual values of the expansion coefficients are again scalar products 
of vectors of numbers in the last three columns of Table 2 with the vector of coefficients 
which forms the first column. Note that the numbers presented in Table 2 correspond to 
the expansion coefficients of (p± , p2 , P3 ; e) /C| . 

7 Conclusions 



This paper finishes the calculation of the integrated doubly unresolved approximate cross 
section of the NNLO subtraction formalism of Refs. [22, 23], and thus completes the defi- 
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nition of the subtraction scheme. 

In particular, here we computed the double soft-type contributions to the integrated 
doubly unresolved approximate cross section of Ref. [22]. The integrated counterterms 
were evaluated in terms of a set of basic double soft integrals, which were calculated using 
Mellin-Barnes techniques. These contributions represented the last missing ingredients 
needed to evaluate f 2 d<r^'^ 2 , as the collinear pieces of the doubly unresolved subtraction 
terms were integrated in Ref. [42]. The final result can be written as the product (in colour 
space) of the Born cross section times the doubly unresolved insertion operator, I^. We 
were able to compute this insertion operator analytically up to 0(e~ 2 ), while the rest of 
the expansion coefficients were evaluated numerically. 

As stressed above, the definition of the NNLO subtraction formalism of Refs. [22, 23] 
is now complete, and the evaluation of the finite doubly virtual cross section of Eq. (1.5) 
is feasible for electron-positron annihilation into two- and three-jets. (A few integrals 
were evaluated specifically for three-jet kinematics, so for a higher number of jets, some 
more work is required.) In particular, all integrated approximate cross sections appearing 
in Eq. (1.5) are known analytically up to 0(e -2 ) and hence, the cancellation of poles in 
d°m NLO ma y be checked explicitly to this order. Indeed, we have checked that the 1/e 4 and 
1/e 3 poles cancel for e + e" — > 2, 3 jets, independently of the value of j/o- The cancellation 
of the subleading poles must be checked numerically, and the details will be presented 
elsewhere. 
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A Modified double soft-type subtraction terms 

In previous publications, we have presented an easy modification to the NNLO subtraction 
scheme of Refs. [22, 23]. The parts of these modifications relevant to all singly unresolved 
approximate cross sections, i.e., dcr^'^ 1 in Eq. (1.3) as well as dcr^'f 1 and ( J x do-J^' 2 1 ) 1 



in Eq. (1.4), were given in Ref. [25], while the parts relevant to the iterated singly unresolved 
approximate cross section, dcr^'^ 12 in Eq. (1.3), were spelled out in Ref. [32]. Finally, the 
modified doubly unresolved approximate cross section was discussed in Ref. [42]. 

The introduction of this modification serves two purposes. First, it makes the inte- 
grated subtraction terms independent of m, the number of hard partons, see Refs. [25, 32, 
42] for a detailed discussion. Second, it allows to constrain the subtractions to near the 
singular regions in phase space, which improves the efficiency of the numerical implemen- 
tation. 

In this paper, we only need the precise definition of the modified double soft-type 
subtraction terms. We recall from Ref. [42] that these are obtained from the original sub- 
traction terms in Ref. [22] via multiplication by the simple factor of f(yo,y r Q + UsQ — 
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Double soft-type counterterms 


Subtraction term 


Momentum mapping 


Function 


o(0,0) rq -(0,0) 

r c(o,o) r <j(o,o) 

n no o(0i0) 
, -w-.s u;, ir;s< ;, »' , s 


Mm+2 ^ 


f(V0, VrQ + VsQ-yrs, d'(m, e)) 



Table 3. The modified double soft-type subtraction terms are obtained from the original coun- 
terterms (first column) by multiplication with f(yo,VrQ + UsQ — Vrsid! (m, e)) (last column). Also 
shown is the momentum mappings used to define the subtraction terms (middle column). 



y rs ,d'(m,e)), where 

f(z ,z,p) = (i-zy p e(z -z). (a.i) 

We also reproduce the relevant part of Table 20 of Ref. [42] here, as Table 3. We note that 
the form of d'(m, e) appearing in Table 3 is actually fixed by the prescription in Ref. [25] and 
the requirement that the modified subtraction terms still correctly regularise all kinematic 
singularities. We have 

d'( m , e) = m(l - e) - d' Q , (A.2) 
where, as in Refs. [25, 32, 42], the parameter d' takes the form 

d' = D' + d[e , (A.3) 

with D' > 2 an integer and d[ real. Throughout this paper, we use the specific choice 
-Dq = 3 and d[ = —3. Finally, we note that in terms of y r Q and y s Q, we have 

f{yo,y r Q + y S Q - y rs , d'(m, e)) = [(l - j/fo)(i - ysQ)] d ' ~ m{1 ~ e) Q{yo - y?Q - y S Q + yfQy s o) , 

(A.4) 

where y rs = y rs and y rQ = (f - y s Q)yfQ + yrs have been used, see Eq. (5.7). 
B Integrating the double soft counterterms 

From Eqs. (4.1)-(4.3) we see that the double soft counterterms involve the integrands 

< SikirjSjiis) , S ik (r, s) , — I 2— j V , (B.l) 

I s rs V s i(rs) s k(rs) S i(rs) ' J 

where Sif-(r, s) is given in Eq. (4.4). 

Let us proceed to identify the independent kinematic structures that we need to inte- 
grate. We begin with the abelian double soft gluon contribution, i.e., the first term in the 
list of Eq. (B.l). Recall that in this term, the indices i, j, k, I are constrained only by the 
requirements that i ^ k and j ^ I, but there is no restriction on whether or not i, k is 
equal to j, I. Hence, we have the following situations: (i) all of i, k, j and I are distinct, (ii) 
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only three of the four indices are distinct, e.g., I = k and (iii) only two indices are distinct, 
e.g., j = i and I = k. 

For case (i) to occur, there must be at least four hard patrons in the process. Obviously 
this does not happen in the calculation of two- and three-jet observables, and we will not 
consider it in this paper. For cases (ii) and (iii), we have simply 



S ik (r)S jk (s) = and S lk (r)S ik (s) = ^ . (B.2) 

SirSkrSjs$ks SirSkr$is$ks 



Turning to the non-abelian part of the double soft gluon formula, the second term in 
the list of Eq. (B.l), we make two observations. First, the factorised phase-space measure, 
[dp^miPriPs) Q)]i is clearly symmetric under the exchange of pr and ps (see Eq. (2.10)). 
Second, the double soft current appears in the subtraction terms summed over its indices, 
i and k in this case. Hence, we are free to exchange i and k in individual terms in Sik(r, s) 
without changing the value of the total integrated subtraction term, since % and k are 
merely summation indices. Using the freedom to make r -H- s and/or % k replacements 
whenever convenient, with some algebra 4 we can derive the following form of 5jfc(r, s), 

Sik(r, s) -> - + (1 



SisSkr^rs 2 Si r SkrSis$ks \S rs ^i(rs)^k(rs)^rs J ^i(rs)^k(rs)^rs 



+ 2- 



^i(rs)^k(rs)^ir^kr $i(rs) ^kr^r 



(B.3) 



whose integral is equal to the integral of the original expression. 

Finally, consider the expression for double soft quark-antiquark emission, the last term 
in Eq. (B.l). Making use of the freedom to exchange the summation indices i and k, after 
some algebra we find 

1 ( Sir^ks Skr^is Sik$rs ^SifSig^ ^ ^ Sir^kr &ik ^ ^ ^ir 



-,2 \ 2/2 22' 

s i(rs) s k(rs) S i(rs) ' s i(rs) s k(rs) s rs s i(rs) s k(rs) s rs S i(rs) Sr s 

(B.4) 

As in Eq. (B.3), the notation — > above means that the two forms lead to the same total 
integrated subtraction term. 

Examining Eqs. (B.2)-(B.4), we identify the following independent kinematic struc- 
tures 5 that we must integrate 



1 SikSjk 1 s 2 ik s ik s ik 



Z Si r SkrSjsSks ^ SirSkr^isSks SisSkr^rs ^i(rs)^kr^rs 

&ik Sir^kr 1 ^jj. 



' ' 2 ' 2 ' 2 

Si(rs)Sk(rs)SirSkr ^i(rs)^k(rs)^rs Si(rs)Sk(rs)Srs ^rs ^i(rs)' 



(B.5) 



4 We found the Mathematica package $Apart [50] useful during these manipulations. 
J The construction of this basic set is not unique, and we make no claim that this set is linearly inde- 
pendent. 
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The specific normalisations were chosen for later convenience. Finally we introduce the 
following notation, 

2 



^25,8(e; Vo,d' 



~S7 



16tt 2 



Q 



2c 



Q 



2c 



16lT 2 



Y2c 



T 2S ,3(Y ll o ,e;yo,d' ) = ( ~^~Q 



T 2 s,4(Y-r^ iQ ,e;y ,d' ) = ( ~^—Q 2e 



16tt 2 



16^ 



Q 



2c 



Q 



2c 



l 2 sj(Y~ i% o ,e;y ,d' ) = ( —Q 



IQtt 2 



l2r 



^25,g(e;yo,4) = ( —g— Q 



16vH 



2c 



l2e 



2 Si r Sk r Sj s Sk s 

x /(yo, 2/rQ + ysQ - yrs, d'(m, e)) , 

(B.6) 



^ ir kr is°ks 

x /(yo, yrQ + y S Q - y rs , d'(m, e)) , 

(B.7) 



1 s 



ik 



x f(V0, VrQ + VsQ ~ Vrs, d! (m, e)) , 

(B.8) 



wflWiUih — — 

■^(rs^fcr^rs 

x /(yo, yrQ + VsQ - Vrs, d'(m, e)) , 



(B.9) 



[d^PpS^]— — — 

*i(rs)*'A;(rs) b ir i, kr 

x /(yo, yr-Q + y S Q - yrs, <Z'(m, e)) , 



(B.10) 



■"'ik 



Si(rs)Sk(rs)Srs 

x /(yo, yrQ + y S Q - yrs, d'(m, e)) , 

(B.ll) 

x /(yo, yrQ + y S Q - y rs , d'(m, e)) , 

(B.12) 

[dpSjtdp^+J-Jr 

6 rs 

x /(yo, yrQ + y S Q - yrs, d'(m, e)) , 

(B.13) 

6 i(rs) 6 rs 

x /(yo, yrQ + y S Q - y rs , d'(m, e)) . 

(B.14) 
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In terms of these basic integrals, the double soft integrated counterterms are expressed as 
in Eqs. (4.7)-(4.10). 



C Integrating the triple collinear — double soft counterterms 

From Eqs. (4.12)-(4.14), we find that the triple collinear - double soft counterterms involve 
the following integrands 



I Zi,rs ( 1 



2^? 1 r / 1 

^%rs - 1 - <-i,rs I 1 



(C.l) 



Note that we have not written the r <B> s part of the non-abelian triple collinear - double 
soft gluon counterterm (see Eq. (4.13)) since by the r o s symmetry of the factorised 
phase-space measure, its integral is equal to the integral of term shown. 

Now, we would like to identify a basic set of kinematic structures, which we must 
integrate. At first glance, it would seem like we ought to use the identity Zi ;rs + z r A S + z S Ar = 
1 to eliminate one of the z's. However it turns out that it is more convenient to simply 
keep the integrand essentially in its original form, as given in Eq. (C.l), for reasons we shall 
explain shortly. Then we find that up to an r o s exchange, there are nine basic structures 
to integrate. These are: 

1 ^i,rs ^s,ir Sir^is ^r,is^s,ir 1 ^i,rs 



Sir$is Zr,isZs,ir ^i(rs)^Vs (^T,is ^s,ir)^ ^i(rs)'^rs {.%r,is ^s,ir) 2 ^i(rs)^rs (%r,is %s,ir) 
1 Zi,rs 1 ^i,rs 1 ^i,rs 1 ^i : rs 



(C.2) 



The reason the above integrands are convenient is as follows: clearly the results of all these 
integrals could only depend on x-r = 2pi ■ Q/Q 2 . However, all nine of the above integrands 
are degree zero homogenous functions of pi. Hence, they are actually independent of x^, 
and thus are free of any kinematic dependence. This nice property is clearly lost if we 
decompose e.g., Z{ rs in the numerators as z% rs — 1 — %r,is — Zs.ir and perform partial 
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fractioning. Then we introduce the notation 



167T 2 2 \ 2 [ (?) u , («) , 1 Z lrs 



m+lJ s- z ■ z ■ 



x f{yo, VrQ + VsQ ~ Vrs, d'(m, e)) , (C.3) 

\ °6 / J2 b i(rs) b rs \ z r,is T Z Sj j r ; 

x /(z/0) + y s Q - 2/rs, e)) , (C.4) 

T / 7/ \ ^ 167T 2 ,,2e"\ /" rj (r) irj W 1 s ir s is z r,is z s,ir 
l2SM e '>y0id ) = l^^Q im-JnliM-'m+lH IT 71 TT 72 

\ "e / J2 i(rs) rs \ z r,is t z s,ir) 

x /(yo, yrQ + - yrs, d'(rn, e)) , (C.5) 



Si(rs)Srs i. z r,is z s,ir~) 

x f(yo, VrQ + y S Q - y rs , d'(m, e)) , (C.6) 

167T 2 2 A 2 f Jr) (s) 1 ^ 



x f(yo, y r Q + y S Q - y rs , d'(m, e)) , (C.7) 
^,iB(6;w),d&) = (^Q 2e ) 2 /[dpS;i][dpSt +1 i^ ^ 



;m+lJ sv z ■ 



Si(rs)Sir z r,is( z r,is ~i~ z s,ir~) 

x f(yo, y T Q + y s q - y rs , d'(m, e)) , (Ci 

Si(rs)Sis z r,is( z r,is ~i~ z s,ir) 

x /(yo, yrQ + y S Q - y™, d'(m, e)) , (C.9) 



x 25a 7( e ;yo,^) = (^Q 2e )^[dpS;i][dpSt +1 ] 1 Zi 



Si r S rs Zr,is 

x /(yo, y r q + y S Q - y rs , d'(m, e)) , (CIO) 

' 1671-2 2 f MaJ s ) 1 1 



x 25 , 18 ( e ;yo,<) = (^Q 2e ) ^[dpS;i][drit +1 ] ; 



x f(yo, y r Q + y S Q - y rs , d'{m, e)) . (C.n) 

The triple collinear - double soft integrated counterterms are expressed with these basic 
integrals as in Eqs. (4.15)-(4.17). 

D Integrating the double collinear — double soft counterterm 

Recall from Eq. (4.18) that the integrated double collinear - double soft counterterm in- 
volves the integrand 

J 1 2 z i,r 1 2zj,s \ (D 1) 

^ Sir z r,i Sjs z s,j J 
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We again keep the integrand in its original form, without making use of the identities 
Zi, r + z T) i = 1 and Zj tS + z s j = 1. The reason is clear: since the expression in Eq. (D.l) is 
degree zero homogenous in both pi and pj, there can be no direct dependence on xj or x ~ . 
However, since the phase-space measure is not exactly factorised in pr and pt (it is exactly 
factorised in pr and ps), the integral still depends on Y-rj q. Introducing the notation 

l 2S MY. liQ ,e;yo,d>) = f^VY / [dp%J [dp^ +1 ] \ 

\ / J2 i>i r bj s Z r ^Z s J 

x f(yo, VrQ + VsQ ~ Vrs, d'(m, e)) , (D.2) 
(note the normalisation) we immediately find Eq. (4.19). 

E Integral of the soft collinear — double soft counterterm 

Form Eq. (4.20) we see that the integrated soft collinear - double soft counterterm involves 
the integrand 

where j = (ir) or I = (ir) is also allowed. (Recall that at the level of kinematics, e.g., j = 
(ir) means that the momentum j entering the eikonal factor is p 1 ^ = p^ + Pr-) 

Using the explicit expression of the eikonal factor, Eq. (2.18), we find that Eq. (E.l) 
evaluates as 



(E.2) 

Sir %r,i Sj s Si s Si r Z T: i S^ r ^ s Si s 

where the first expression corresponds to the case when both j and I are distinct from (ir), 
while the second corresponds to j = (ir). We introduce the notation 

1671-2 r> 2e\ 2 fn (r),M (s) i 1 Zi,r Sjl 



Sir Zr,i Sjs$ls 

x f(yo, VrQ + VsQ ~ Vrs, d'(m, e)) , (E.3) 

l^n&V fuJ*)VAJ') 1 1 S (ir)l 



x f(V0, VrQ + VsQ ~ Vrs, d'{m, e)) , (E.4) 



Sir %r,i S(j r ) s Sjs 



and obtain the soft collinear - double soft integrated counterterm as in Eqs. (4.21) and (4.22), 
for j, I ^ (ir) and j = (ir) respectively. 

F Integrating the triple collinear — soft collinear — double soft countert- 
erm 

Recall from Eq. (4.23) that the integrated triple collinear - soft collinear - double soft 
counterterm involves the integrand 

^Zi,rs(Zi,rs ~\~ Z T) is) 



SirS(ir)s%r,isZs,ir 



(F.l) 
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In this case, we will simply split the sum in the numerator and use as our independent 
integrals the following: 



1 



1 Zi 



2 2' 



(F.2) 



With the notation 

^2S,22<X7, e;y ,d' ) 



x f(yo, VrQ + VsQ - Vrs, d! (m, e)) , 
16tt 2 ^\ 2 f r , (f) ir ,» 



S(ir)s s ir Z r,is Z s,ir 



1 s 



i,rs 



1 2 ^(ir)s^ir %s,ir 

x f(yo, y r Q + y S Q - y rs , d'(m, e)) , 



(F.3) 



(F.4) 



(note the normalisation) we immediately find Eq. (4.24). 

G Collinear-type doubly unresolved counterterms to 0(e~ 2 ) 

In this appendix we recall from Ref. [42] the analytic formulae for the flavour-summed 
collinear-type functions in Eq. (3.3). We have 

1. Triple collinear: 



(G.l) 



2. Triple collinear - soft collinear: 



C irs CS ir . s j \x, L/ 



1 2 



In a;, + Z(yo,D' - 1] 



+ o( e - 2 ). 



(G.2) 



3. Double collinear: 



,(o) 

J ir;js 



fifj 



1 



1 



2e 4 2e 3 

+ 0(6- 2 ) . 



2(lnxj + Inxj) - 7/ 4 («f) - 7/,(™f) 



(G.3) 



4. Double collinear - soft collinear: 

+ 0(e~ 2 ) . (G.4) 
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5. Soft collinear: 



GS 



(o) 



12/ \ 1 

i,Y jl;Q ) = + ^ ( In Sj + Z(y ,D' - 1)1 + -3 ( In Y^q - 7 /i (n f )) 

+ 0(e" 2 ), (G.5) 



for j 7^ i, and 



1 2 , 

^4 + ^3 ( Ins, + T,(y ,D - 1 



+ i(lny« ;0 - | 7 /,(n f )) + 0(6" 2 ) . (G.6) 



for j = i. 
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